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1.0  INTRODUCTION 

The  purpose  of  the  heat  transfer  data  base  system  is  to  organize  the  massive 
data  base  generated  under  Contract  F33615-77-C-203Q  and  to  provide  a  conven¬ 
ient  way  of  accessing  the  data  through  the  computer.  The  data  base  system 
solves  the  problem  of  storing  large  volumes  of  printed  data,  since  data  are 
stored  on  magnetic  tape.  The  data  base  is  read  from  tape  into  a  direct-access 
disk  file  that  may  be  accessed  in  a  number  of  ways  using  the  software  routines 
provided.  The  direct  computerized  access  avoids  the  human  errors  associated 
with  transferring  large  quantities  of  data  from  printed  sheets  to  punched 
cards  or  other  computer  input  formats. 

Figure  1  is  a  schematic  of  the  endwall  heat  transfer  data  base  system  at  De¬ 
troit  Diesel  Allison.  The  figure  highlights  the  parts  of  the  data  base  de¬ 
livered  and  described  in  this  volume. 
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Figure  1.  Endwall  heat  transfer  data  base  system 


2.0  SYSTEM  OVERVIEW 


2.1  DATA  CONTENT 

The  concents  of  Che  data  base  include  all  raw  and  reduced  data  and  geometry 
information  used  to  generate  the  results  reported  in  Volume  II  of  this  report. 
All  temperatures  are  given  in  °F,  all  pressures  in  psia,  and  all  dimensions 
(including  relative  coordinate  locations)  in  inches.  Velocities  are  given  in 
ft/sec,  and  Reynolds  numbers  are  based  on  true  chord.  Any  other  unique  di¬ 
mensions  are  stored  in  the  list  program  and  may  be  obtained  from  an  index 
listing. 

The  listing  program  specifies  the  basis  of  several  derived  quantities  in  terms 
of  subscripted  variable  names.  Table  I  shows  the  subscript  convention  and  de¬ 
fines  the  subscripted  variable  names.  Definitions  for  the  derived  loss  para¬ 
meters  also  are  included  in  Table  I. 


TABLE  I.  SUBSCRIPTED  VARIABLE  DEFINITIONS 


Subscript  convention:  1 — measured  inlet 

2 —  measured  exit 
X — ideal  exit 

3 —  uniform  mixed-out  exit 


Variable 


Definition 


PT1 

PS1 

TT1 

PTX 

PSX 

TTX 

PT2 

PS2 

TT2 

0MEGA2 

DS0CP2 


EBAR2 


PT3 

PS3 

0MEGA3 

DSCP3 


EBAR3 


Measured  inlet  total  pressure 

Average  of  measured  inlet  sidewall  static  pressures 
Measured  inlet  total  temperature 
Ideal  exit  total  pressure  (*  PT1) 

Average  of  measured  exit  sidewall  static  pressures 
Ideal  exit  total  temperature  (  =  TT1) 

Exit  cone  probe  total  pressure 
Exit  cone  probe  static  pressure 
Exit  probe  total  temperature 
(PT1-PT2)/(PTX-PSX) 

Ln  (TT2/TTl)-^y^)  Ln  (PT2/PT1 ) 


Zzl 

(PTX/PT2)  y  -1. 
— 

(PTX/PS2)  7  -1. 


Calculated  mixed-out  total  pressure 
Calculated  mixed-out  static  pressure 
(PTX— PT3)/(PTX— PS3) 


Ln  (TTX/TT1) 


Ln  (PT3/PT1) 


Zzl 

(PTX/PT3)  7  -1. 


(PTX/PS3)  7  -1. 
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Most  of  the  raw  data  Items  stored  In  the  data  base  have  been  reviewed  and 
edited  to  flag  obvious  bad  readings  due  to  instrumentation  failures.  The 
edited  readings  were  flagged  by  replacing  the  bad  data  with  the  value  -1*  All 
thermocouple  readings  for  the  heat  transfer  endwall  have  been  edited  to  in¬ 
clude  only  those  readings  used  as  input  to  the  finite  element  analysis  pro¬ 
gram.  Thermocouple  data  from  the  adiabatic  endwall,  the  vane,  and  upstream 
total  temperature  rakes  were  also  edited.  Upstream  total  and  static  pres¬ 
sures,  endwall  static  pressures,  and  exit  static  pressures  were  edited  as  well. 

Data  associated  with  the  boundary  layer  pressure  and  temperature  rakes  were 
not  edited.  The  data  from  these  rakes  were  reliable  for  early  runs,  but  indi¬ 
cated  significant  instrumentation  failures  as  time  went  on.  Since  data  from 
the  traversing  inlet  probe  have  been  supplied  to  provide  boundary  layer  infor¬ 
mation,  the  boundary  layer  rake  data  were  not  edited.  These  data  have  been 
included  in  the  data  base,  however,  and  should  therefore  be  used  with  caution. 

2.2  DAT A  0  RGANI ZAT 10  N 

The  data  base  consists  of  sets  of  data,  and  each  set  includes  all  data  with  a 
common  identifier.  The  present  data  base  consists  of  two  types  of  data  sets, 
but  may  be  expanded  as  the  user  desires. 

The  first  type  of  data  set  in  the  data  base  is  the  single  data  set  GEOMETRY. 
This  set  includes  relative  coordinates  for  all  endwall  instrumentation  and 
vane  thermocouple s.  The  data  set  also  includes  a  vane  profile  description. 

Each  item  included  in  the  set  GEOMETRY  has  an  assigned  identification  number 
that  is  used  in  accessing  the  data.  Table  II  lists  each  GEOMETRY  variable  and 
its  identification  (I.D.)  number. 


TABLE  ti.  DATA  SET  GEOMETRY  VARIABLE  l.D.  MAP 


ID 

NO. 

DESCRIPTION 

DIMENSIONS 

PRIMAR. 

SECONDAR' 

1 

H.T.  ENDWALL  HOT  SIDE  T/C  COORDINATES 

2 

53 

2 

H.T.  ENDWALL  COLD  SIDE  T/C  COORDINATES 

2 

31 

3 

FINITE  ELEMENT  MODEL  NODE  COORDINATES 

2 

338 

4 

ADIABATIC  ENDWALL  T/C  COORDINATES 

2 

90 

5 

ENDWALL  STATIC  PRESSURE  COORDINATES 

2 

42 

6 

VANE  T/C  COORDINATES 

2 

113 

7 

ENDWALL  BOUNDARY  POINTS 

3 

85 

8 

VANE  PROFILE  COORDINATES 

2 

60 

9 

VANE  TRUE  CHORD 

1 

1 

10 

VANE  SPACING 

1 

1 

11 

VANE  LEADING  EDGE  RADIUS 

1 

1 

12 

VANE  TRAILING  EDGE  RADIUS 

1 

1 

13 

VANE  SETTING  ANGLE 

1 

1 

14 

CASCADE  THROAT  HEIGHT 

1 

1 

10 
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The  second  type  of  data  set  currently  stored  in  the  data  base  is  made  up  of  all 
data  associated  with  a  particular  experimental  run  number*  Forty-four  such 
data  sets,  RUNXXXXX,  are  stored  in  the  data  base-  Each  variable  stored  for  a 
RUNXXXXX  set  has  an  associated  identification  number.  The  general  map  of  var¬ 
iables  and  associated  identification  numbers  used  for  all  RUNXXXXX  data  sets 
is  shown  in  Table  III.  See  Section  3  of  Volume  II  of  this  report  to  determine 
the  run  conditions  and  data  available  for  a  given  run  number.  A  detailed  in¬ 
dex  of  the  variables  stored  for  each  run  number  may  also  be  generated  using 
the  program  BLIST. 

2.3  MANAGEMENT  SOFTWARE 

The  management  software  provided  with  the  data  base  includes  two  stand-alone 
programs  for  editing  and  listing  data,  two  interpolation  function  subroutines, 

|  and  18  general  access  subroutines  for  data  item  retrieval  and  storage. 

The  stand-alone  programs  are  MAINT  and  BLIST.  These  programs  use  the  access 
subroutines  to  perform  routine  tasks.  MAINT  facilitates  the  interrogation  of 
the  data  base  with  the  option  of  editing,  creating,  or  deleting  data  sets. 

BLIST  provides  for  listing  all  or  any  subset  of  the  data  base. 

The  interpolation  functions  have  been  provided  for  the  analyst  to  facilitate 
use  of  the  data  base.  ANTRP  provides  a  single  value  for  any  variable  asso¬ 
ciated  with  the  endwall  region  for  a  given  (x,y) .  BNTRP  returns  a  single  val¬ 
ue  of  any  variable  associated  with  the  traversing  exit  probe  for  a  given  (z,y) 
point  in  the  defined  exit  plane.  The  interpolation  functions  may  be  incor¬ 
porated  in  any  user-written  program  accessing  the  data  base. 

The  file  access  subroutines  are  the  building  blocks  for  programs  interacting 
with  the  data  base.  They  are  generalized  subroutines  for  item  storage  and 
retrieval.  The  access  subroutines  are  written  in  standard  FORTRAN,  and  may  be 
incorporated  in  any  program  requiring  data  base  interaction. 
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TABLE  III.  DATA  SET  RUNXXXXX  VARIABLE  I.D 
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3.0  USER'S  MANUAL 


3.1  SOFTWARE  DESCRIPTIONS/USAGE  INSTRUCTIONS 

3.1.1  Main  Programs 

The  main  programs  supplied  with  the  data  base  are  a  maintenance  program  and  a 
bulk  data  listing  program.  Both  are  written  to  read  the  data  base  from  mag¬ 
netic  tape  and  create  a  random-access  data  bank  disk  file.  Both  programs  read 
input  data  from  punched  cards  and  output  information  to  a  line  printer  in  133- 
character  records.  The  general  input  deck  structure  is  outlined  in  Section 
3.2.  The  programs  interact  with  the  random-access  data  bank,  processing  a 
series  of  commands.  The  maintenance  program  creates  a  new  data  base  magnetic 
tape  file. 

Program  1 — Maintenance  (MAINT) 

The  purpose  of  the  maintenance  program,  MAINT,  is  to  provide  a  means  of  making 
minor  changes  in  the  data  base.  It  is  not  intended  to  be  used  for  entering 
large  volumes  of  data.  The  program  generates  a  printed  record  of  changes  made 
as  an  audit  trail,  listing  old  and  new  values.  The  program  also  permits  in¬ 
quiry  about  the  value  of  any  data  item  without  updating  it. 

The  maintenance  program  begins  by  reading  a  sequential  tape  file  of  the  data 
base  and  creating  a  random-access  data  bank  disk  file.  The  program  then  com¬ 
municates  only  with  the  random-access  data  bank,  processing  a  series  of  input 
commands  sequentially.  When  the  command  list  is  exhausted,  the  data  bank  disk 
file  is  closed  and  copied  to  a  new  sequential  tape  file,  forming  the  new  data 
base. 

Input  to  MAINT  consists  of  a  series  of  commands  using  key  words  that  may  be 
arranged  to  suit  the  user's  needs.  The  format  of  the  commands  is  shown  in 
Figure  2.  Each  command  line  is  30  characters  long.  Key  words  are  underlined, 
and  when  shown  as  beginning  in  column  one,  they  must  appear  there.  When 
blanks  are  underlined,  they  are  an  integral  part  of  the  key  word. 

All  numbers  are  read  in  free  format.  Decimal  points  and  exponential  notations 
may  be  used,  as  in  the  FORTRAN  language.  Repeated  values  may  be  coded  as  N  * 
#.//#,  where  N  is  the  number  of  occurrences.  All  free-format  numbers  must  be 
separated  by  blanks  or  commas.  Data  may  be  continued  on  subsequent  lines,  and 
may  begin  on  the  line  following  the  identification  specification.  No  number 
field  may  be  broken  by  a  continuation. 

The  commands  shown  in  Figure  2  can  be  separated  into  four  major  types.  The 
mode  command  specifies  the  function  to  be  performed,  the  set  command  specifies 
the  data  set  to  be  accessed,  the  l.D.  commands  specify  the  data  manipulation 
within  the  data  set,  and  the  end  command  signals  the  end  of  the  maintenance 
session. 
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Figure  2.  Maintenance  program  commands. 
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The  mode  command  specifies  the  function  to  be  performed  by  subsequent  com¬ 
mands.  In  the  UPDATE  mode,  the  data  base  is  changed  as  dictated  by  following 
commands.  The  program  begins  in  the  UPDATE  mode.  The  INQUIRY  mode  will  re¬ 
turn  the  data  stored  under  the  identification  following. 

The  mode  may  be  changed  at  any  point  in  the  run,  and  the  program  remains  in 
that  mode  until  another  mode  command  is  encountered.  In  the  INQUIRY  mode,  the 
program  ignores  any  data  encountered. 

The  second  major  type  of  command  is  used  to  specify  the  data  set  under  which 
the  following  identification  numbers  are  to  be  found  or  inserted.  The  com¬ 
mands  are  numbered  3  and  4  in  Figure  2.  Data  stored  in  the  data  base  are 
divided  into  data  sets.  All  data  for  a  given  experimental  run  number  are 
stored  under  a  data  set  called  RUNXXXXX,  where  the  Xs  represent  the  run  num¬ 
ber.  Instrumentation  coordinates  and  other  geometry  information  are  stored  in 
a  set  called  GEOMETRY.  The  SET  command  is  used  to  enter  the  eight -character 
set  name  directly.  The  RUN  command  synthesizes  the  eight-character  set  name 
for  the  free-format  run  number  given,  and  then  performs  the  same  action  as  the 
SET  command.  One  of  the  option  key  words  may  be  included  (beyond  column  12  on 
the  command  line)  with  either  of  these  commands.  DELETE  causes  the  entire 
set,  including  all  of  its  data,  to  be  deleted  from  the  data  base.  CREATE  must 
be  present  on  the  set  command  line  (beyond  column  12)  if  a  new  data  set  is  to 
be  initiated.  The  default  action  of  the  SET  command  is  designed  to  use  a  set 
if  it  exists,  and  to  print  an  error  message  if  it  does  not,  unless  CREATE  is 
specified . 

The  third  major  type  of  command,  the  I.D.  command,  is  shown  on  lines  5  through 
14  of  Figure  2.  The  I.D.  command  is  used  to  update  or  inquire  about  the  value 
of  an  individual  data  item  within  a  data  set,  depending  on  the  most  recent 
mode  specification.  The  item  specified  may  be  a  scalar,  an  entire  vector,  an 
element  of  a  vector,  an  entire  matrix,  an  element  of  a  matrix,  a  row  of  a 
matrix,  or  a  column  of  a  matrix. 

When  the  program  is  in  the  UPDATE  mode,  the  list  of  new  data  is  entered  ac¬ 
cording  to  the  identification  number  and  subscript  specification,  as  shown. 
Since  the  list  may  be  continued  on  subsequent  lines,  the  program  will  not  con¬ 
sider  the  command  to  be  complete  until  it  sees  the  beginning  of  the  next  com¬ 
mand.  If  the  program  is  running  in  the  INQUIRY  mode,  any  data  list  erroneous¬ 
ly  supplied  will  be  ignored. 

Commands  numbered  5  through  11  in  Figure  1  illustrate  the  various  data  manipu¬ 
lations  possible  using  MAINT.  Command  5  shows  how  a  scalar  would  be  entered 
in  the  UPDATE  mode.  In  the  INQUIRY  mode,  no  data  should  be  present  on  the 
command  line.  Command  6  shows  how  to  refer  to  an  element  of  a  vector,  with 
the  subscript  placed  in  parentheses,  as  shown.  This  command  is  normally  used 
for  vectors  that  already  exist  in  the  data  base.  If  a  new  vector  is  being 
created,  however,  it  will  be  initialized  with  maximum  dimensions  equal  to  the 
subscript,  and  the  last  element  will  be  assigned  the  data  supplied. 

Command  7  shows  how  to  update  or  interrogate  an  entire  vector.  The  paren¬ 
theses  and  asterisk  are  not  really  required  for  this  case.  If  a  new  vector  is 
being  created,  the  list  length  determines  the  maximum  dimensions  it  is  as¬ 
signed  in  the  data  base.  When  an  existing  vector  is  overstored,  the  length  of 
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the  list  may  not  exceed  the  dimensions  previously  assigned  to  that  vector.  To 
increase  the  dimensions  of  the  vector,  the  vector  must  first  be  deleted  and 
then  re-created.  When  the  list  length  is  less  than  the  length  of  the  existing 
vector,  excess  elements  are  not  updated. 

Commands  8  through  11  show  the  syntax  for  updating  an  element,  a  column,  a 
row,  or  an  entire  matrix.  The  use  of  command  11  requires  that  the  matrix  must 
have  been  previously  initialized,  so  that  its  primary  dimensions  may  be  ob¬ 
tained  from  the  data  base.  Command  8  is  used  to  initialize  a  new  matrix  with 
the  maximum  dimensions  given.  As  in  the  case  of  a  vector,  the  maximum  dimen¬ 
sions  may  not  be  exceeded.  The  matrix  must  be  deleted  and  re-created  to  ex¬ 
pand  its  size.  Commands  9  and  10  show  how  to  update  a  row  and  column  of  a 
matrix,  respectively.  If  the  data  list  is  shorter  than  the  matrix  row  or  col¬ 
umn  previously  initialized,  the  update  starts  at  the  beginning  of  the  row  or 
column  and  proceeds  until  the  list  is  exhausted.  The  remaining  elements  are 
not  updated.  When  commands  8  through  11  are  used  in  the  INQUIRY  mode,  any 
data  list  is  ignc.Jd. 

Command  12  in  Figure  1  illustrates  a  data  base  feature  that  facilitates  entry 
of  tabular  data.  The  I.D.s  specify  the  headings  of  the  table  and  the  data  are 
input  as  ID1  (1),  ID2  (1),  ID3  (1),  ID1  (2),  ID2  (2),  and  so  on.  The  dimen¬ 
sions  must  match  for  all  I.D.s,  and  the  length  of  the  data  list  must  be  a  mul¬ 
tiple  of  the  number  of  I.D.s  specified.  The  subscript  information  specified 
in  parentheses  is  not  required  for  vectors,  but  must  be  supplied  for  matrices. 

Commands  13  and  14  in  Figure  1  show  how  to  delete  a  variable  or  variables. 

Once  dimensions  have  been  set  for  a  variable,  they  are  permanently  associated 
with  the  I.D.  and  may  not  be  changed  unless  the  I.D.  is  first  deleted  by  means 
of  one  of  these  commands. 

The  final  command  Is  numbered  15  in  Figure  1.  The  END  command  is  used  at  the 
end  of  the  command  string  to  trigger  the  closing  of  the  files  and  data  base 
and  the  orderly  termination  of  the  maintenance  session. 

A  sample  input  data  set  for  MAINT  is  shown  in  Table  IV.  The  output  print  gen¬ 
erated  for  this  sample  maintenance  session  is  shown  in  Tabic  V. 

The  sample  session  illustrates  several  important  rules  for  the  maintenance 
program.  In  command  3,  note  that  I.D.  1  must  be  deleted  before  its  maximum 
dimensions  can  be  altered.  Command  4  shows  that  a  vector  size  need  not  be 
specified  for  initialization.  Command  7  shows  the  input  convention  for  ele¬ 
ments  of  a  matrix*  The  primary  index  varies  fastest.  Commands  11,  12,  and  13 
illustrate  the  flexibility  of  free-format  input. 

In  the  inquiry  mode,  note  the  action  resulting  from  commands  17  and  18.  Com¬ 
mand  17  calls  out  a  matrix  I.D.,  but  only  one  column  of  the  matrix  is  returned 
if  (*,*)  is  not  present  on  the  command  line.  Commands  20  through  26  recall 
the  data  stored  above  and  23  illustrates  the  output  for  multiple  I.D.  lists 
supplied  in  the  INQUIRY  mode.  Commands  27  through  31  demonstrate  the  use  of 
the  DELETE  command  and  the  error  message  printed  when  a  deleted  variable  is 
called. 
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TABLE  IV.  SAMPLE  INPUT  FOR  MAINTENANCE  PROGRAM 


SC  T  E 
IP  1 
ID  1 
ID  I 
ID  1 
ID  2 
ID  2 
ID  3 
ID 

ID  5 
ID  3 
3  4  « 
3  4  1 

inou: 
run  : 
ID  1 
ID  2 
ID  2 f 
ID  25 
SCT  t 
ID  1 
in  2 

ID  3 
ID  3 
ID  4 
ID  ? 
UPO&l 
ID  4 
inqu: 
ID  4 
ID  5 
END 


(AMPLE 
193  1 

>  2^.  : 
:  *  j  7  i 
;  2 » 3 1  : 

[4)  1  . 
4)  l  . 
14  )  1  . 

ID  4  I  ■ 


DELETE 
1731 
4  17  7>) 

1 . 

11  21  12 


5  (*) 


1  3  23 


5 

5  <*,*) 
EXAMPLE 


delete 


Program  2 — Bulk  Data  Listing  (BLIST) 

The  purpose  of  the  listing  program,  BLIST,  is  to  provide  a  full  or  partial 
list  of  data  base  contents.  The  program  provides  for  the  listing  of  entire 
arrays  or  simply  an  index  of  variables  stored  for  a  given  data  set. 

The  BLIST  software  is  comprised  of  a  mainline,  several  "title"  subroutines, 
and  a  generalized  listing  subroutine. 

The  mainline  program  first  generates  the  random-access  data  bank  from  a  se¬ 
quential  input  tape  file  of  the  data  base.  It  then  proceeds  to  read  and  in¬ 
terpret  the  input  data  commands.  The  mainline  processes  the  input  commands 
sequentially  and  calls  subroutines  to  generate  the  titles  and  variable  lists 
specified.  For  each  variable  list,  the  code  consists  of  a  title  subroutine 
call,  followed  by  a  listing  subroutine  call.  When  the  END  command  is  en¬ 
countered,  the  mainline  closes  the  random-access  data  bank.  BLIST  does  not 
generate  a  new  sequential  tape  file. 

The  "title"  subroutines  include  subroutines  GTITL  and  TITL1  through  TITL10. 
GT1TL  retrieves  the  descriptive  variable  name  corresponding  to  a  specified 
GEOMETRY  I.D.  number.  TITL1  through  T1TL10  contain  the  descriptive  variable 
names  for  I.D.  numbers  of  data  sets  of  the  type  RUNXXXXX.  If  new  variable 
I.D.s  are  added  to  the  data  base,  additional  title  subroutines  must  be  gen¬ 
erated  for  BLIST. 


19 


TABLE  V.  SAMPLE  OUTPUT  FOR  MAINTENANCE  PROGRAM 


OBANK  LO  AC  COMPLETED 
#**SET  EXAMPLE  CREATE 
***ID  1  1981 

DATA  SET  EXAMPLE  ID  NUMBER  i  SIZE 
ELEMENT  OLD  VALUE  NEW  VALUE 

Uli'l  *’  -“delete 

***ID  1  5  26 »  1961 

DATA  SET  EXAMPLE  ID  NUMBER  1  SIZE 
ELEMENT  OLD  VALUE  NEW  VALUE 

(  1,  1)  .0  5.00000 

[  3,  lj  .0  1931.00 

*4*ID  1  <*)  7  A  1776 


DATA  SET  EXAMPLE  ID  NUMBER  1  SIZE 
ELEMENT  OLD  VAi UE  NEW  VALUE 

(  1,  1)  5.00000  7.00000 

(  3,  1)  1981.00  1776.00 

***I0  2  (2,3)  1. 

DATA  SET  EXAMPLE  ID  NUMBER  2  SIZE 
ELEMENT  OLD  VALUE  NEW  VALUE 

<  2,  3)  .0  1.00000 

***I0  2  (*,*)  11  21  12  22  13  23 

DATA  SET  EXAMPLE  ID  NUMBER  2  SIZE 
ELEMENT  OLD  VALUE  NEW  VALUE 

(  1,  1)  .0  11.0000 

(  1,  2)  1981,00  12.0000 

(  1,  3)  .001581E-78  13.0000 

*4*10  3  (4)  1. 

DATA  SET  EXAMPLE  ID  NUMBER  3  SIZE 
ELEMENT  OLO  VALUE  NEW  VALUE 

(  4,  1)  .0  1.00000 

* ** 1 0  4  (4)  1. 

DATA  SET  EXAMPLE  ID  NUMBER  4  SIZE 
ELEMENT  OLD  VALUE  NEW  VALUE 

(  4,  1)  .0  1.00000 

*4*10  5  (4)  1. 

CATA  SET  EXAMPLE  ID  NUMBER  5  SIZE 
ELEMENT  OLO  VALUE  NEW  VALUE 

(  4,  1)  .0  1.00000 

***ID  3  ID  4  10  5  (*)  3  4  5  3 

***3  4  5 

***3  4  5 

DATA  SET  EXAMPLE  ID  NUMBER  3  SIZE 
ELEMENT  CLO  VALUE  NEW  VALUE 

(  1.  1)  .0  3.00000 

(  3,  1)  1981.00  3.00000 

OATA  SET  EXAMPLE  ID  NUMBER  4  SIZE 
ELEMENT  OLO  VALUE  NEW  VALUE 

(  1.  1)  .0  4.00000 

1  3,  1)  1981.00  4.00000 


(  3,  1) 

ELEMENT  OLC  VALUE 

(  2,  1)  .0 


(  3,  1) 

ELEMENT  OLD  VALUE 
(  2,  I)  2O.000G 


NEW  VALUE 

26.0000 


NEW  VALUE 
4.00000 


<  2,  3) 


2,  3) 

ELEMENT  OLC  VALUE 
2,  1)  26.000C 

2,  2  .404820E- 

2,  3)  1.00000 


NEW  VALUE 
21 .0000 
78  22.0000 

23.0000 


(  4,  1) 


4,  l) 

ELEMENT  OLD  VALUE 
(  2,  1)  26.0000 

(  4,  1)  1.00000 

4,  1) 

ELEMENT  OLD  VALUE 
(  2,  1)  26.0000 

(  4,  1)  1.00000 


NEW  VALUE 
3.00000 
3.00000 


NEW  VALUE 
4.00000 
4.00000 


DATA  SET  EXAMPLE  10  NUMBER  5  SIZE  ( 
ELEMENT  DEVALUE  NEWVALU^ 

(  3.  1)  1981.00  5.00000 

***in3uiry 

•**RUN  122 
***ID  1 

DATA  SET  RUN  122  ID  NUMBER  1  SIZE  ( 
ELEMENT  OLO  VALUE  NEW  VALUE 
<  1,  1)  1.00000 

1)  80.0000 


***l3’i 


l )  OL^lfobHo 
1)  1.00000 


eCeMEn!  OLO  VALUE 
(  2,  1)  7.00000 


NEW  VA 
5.00 
5.00 


. 00000 


NEW  VALUE 


TABLE  V.  (GONT) 


OATA  SET  RUN  122  ID  NUMBER  2  SIZE  < 
ELEMENT  OLD  VALUE  NEW  VALUE 

<  1.  1)  769.180 

***10  25 

DATA  SET  RUN  1 22  ID  NUMBER 
ELEMENT  _  OLO  VALUE  Nl 


1) 


»  1, 

|  3. 

(  5, 

***ID  25 


l! 

li. 


*) 


59.3718 

48.2730 

59.9812 


ATA  S 

El|M 

**SEf 
**I  D 


ET  RUN 
ENT 

1 
1 

l 
I! 

EXAMPLE 

1 


122  ID  NUMBER 
OLD  VALUE  Nf 
59.3716 
48.2730 
59.9812 
59.8329 
59.8146 
59.9251 


DATA  SET  EXAMPLE  ID  NUMBER 
ELEMENT  OLD  VALUE  NEW 

(  1,  1)  7.00000 

(  3,  1)  1776.00 

***ID  2 

DATA  SET  EXAMPLE  ID  NUMB 67 
ELEMENT  OLO  VALUE  NEW 

(  I,  1)  11.0000 

*«*ID  3  ID  4  ID  5 

DATA  SET  EXAMPLE  ID  NOMD/.fi 
ELEMENT  OLO  VALUE  *'EW 

(  1,  1)  3.00000 

(  3,  1)  3.C000P! 

DATA  9ET  EXAMPLE  ID  NUM8ER 
ELEMENT. .  CLO  VALUE  NEW 


I. 

3, 


il 


4.00000 

4.00000 


DATA  SET  EXAMPLE  ID  NUMBER 
ELEMENT  .  OLD  VALUE 


««*I 


1. 

8*: 


i! 


5.00000 

5.00000 


DATA  SET  EXAMPLE  ID  NUMBER 
ELEMENT  CLD  VALUE 


3.00000 

3.00000 


Lj;,  i! 


DATA  SET  EXAMPLE  ID  NUMBER 
ELEMENT  OLD  VALUE 


i* 

>**ID5 


1) 


4.00000 

4.00000 


DATA  SET  EXAMPLE  10  NUMBER 
ELEMENT  OLD  VALUE  N 
(  It  1)  5.00000 

1  3 1  1)  5.00000 

***ID  4^0  5  DELETE 
♦♦♦INQUIRY 
***I0  4 

OUNABLE  TO  LOCATE 
CCMMANQ  ERROR 
***ID 


25  SIZE 

< 

5.  2) 

VALUE 

ELEMENT 

OLD  VALUE 

NEW 

VALUE 

l 

2, 

1) 

59.8754 

( 

4, 

n 

59.9444 

25  SIZE 

( 

5.  2) 

VALUE 

ELEMENT 

OLO  VALUE 

NEW 

VALUE 

( 

2  , 

1) 

59.8754 

( 

4, 

1) 

59.9444 

( 

2, 

2) 

59.8120 

1 

4, 

2) 

53.8518 

1  SIZE 
VALUE 

( 

Ii!emenI 

OLD  VALUE 

NEW 

VALUE 

( 

2. 

1) 

4.00000 

2  SIZE 

( 

2.  3) 

VALUE 

ELEMENT 

OLD  VALUE 

NEW 

VALUE 

( 

2. 

1) 

21.0000 

3  SIZE 

( 

4,  1) 

value 

ELEMENT 

OLD  VALUE 

NEW 

VALUE 

( 

2. 

1) 

3.00000 

( 

4, 

ll 

3.00000 

4  SIZE 

{ 

4.  1) 

VALUE 

ELEMENT 

OLD  VALUE 

NEW 

VALUE 

{ 

2, 

1) 

4.00000 

( 

4,' 

ll 

4.00000 

5  SIZE 

( 

4,  1) 

VALUE 

ELEMENT 

OLD  VALUE 

NEW 

VALUE 

( 

2, 

1) 

5 .00000 

( 

4, 

li 

5.00000 

3  SIZE 
VALUE 

( 

elemenI 

OLD  VALUE 

NEW 

VALUE 

< 

2  . 

1) 

3.00000 

( 

4, 

1) 

3.00000 

4  SIZE 
VALUE 

< 

4,  1) 

ELEMENT 

OLD  VALUE 

NEW 

VALUE 

( 

2  . 

1) 

4.00000 

( 

4, 

li 

4.00000 

5  SIZE 
VALUE 

( 

4.  1) 

ELEMENT 

OLO  VALUE 

NEW 

VALUE 

( 

2. 

1) 

5.00000 

( 

4, 

1) 

5.00000 

4  IN  SET  EXAMPLE 


OUNABLE' TO  LOCATE 
COMMAND  ERROR 

***ENO 

END  OF  UPOATES 


5  IN  SET  EXAMPLE 


The  general  listing  subroutine,  WRIT,  prints  the  variable  description  obtained 
from  the  title  subroutine  and  then  retrieves  and  prints  the  data  arrays  if  the 
program  is  in  the  ARRAY  mode*  If  the  program  is  in  the  INDEX  mode,  only  the 
title  line  is  printed.  Subroutine  WRIT  is  general  and  may  be  used  to  list  new 
variables  added  to  the  data  base,  providing  the  title  Information  is  correctly 
supplied . 

Input  to  BLIST  consists  of  a  series  of  commands  using  key  words.  The  commands 
may  be  arranged  to  list  all  or  any  subset  of  the  data  base.  The  commands  for 
BLIST  are  shown  in  Figure  3.  Each  command  line  is  80  characters  long*  Key 
words  are  underlined,  and  when  shown  as  beginning  in  column  one,  they  must 
appear  there.  When  blanks  are  underlined,  they  are  an  integral  part  of  the 
key  word. 

The  commands  shown  in  Figure  5  are  of  three  major  types.  The  mode  command 
specifies  the  function  to  be  performed  by  the  program,  and  the  set  command 
specifies  the  data  set  or  sets  to  be  listed  with  optional  I.D.  specification. 
The  end  command  signals  the  end  of  the  run. 


The  mode  command,  INDEX,  in  Figure  3  will  cause  the  program  to  list  an  index 
of  variable  names  stored  for  the  following  data  set  names.  In  the  ARRAY  mode, 
variable  names  and  their  associated  data  arrays  are  listed.  The  program  be¬ 
gins  in  the  ARRAY  mode  and  remains  in  that  mode  until  another  mode  command  is 
encountered.  The  program  mode  may  be  changed  at  any  point  in  a  listing  series 
of  commands. 

The  second  major  type  of  command  is  the  set  specification  command,  shown  on 
lines  3,  4,  and  5  of  Figure  3.  The  ALL  command  will  cause  all  data  sets 
RUNXXXXX,  to  be  listed  according  to  the  last  mode  specified.  The  RUN  command 
will  produce  a  list  of  all  data  associated  with  the  specified  run  number. 
Geometry  data,  or  an  index  of  geometry  data,  will  be  listed  if  command  5  is 
present  in  the  series  of  commands.  Shown  in  brackets  in  commands  4  and  5  is 
the  optional  I.D.  specification.  If  a  list  of  one  data  array  is  desired,  add¬ 
ing  the  I.D.  number  to  the  set  specification  command  will  produce  the  single 
array  list.  The  I.D.  option  provides  the  flexibility  for  listing  any  subset 
of  the  data  base. 

Finally,  the  END  command  is  used  to  terminate  the  listing  session.  The  com¬ 
mand  triggers  the  orderly  closing  of  files. 

A  sample  input  data  set  for  BUST  is  shown  in  Table  VI.  The  resulting  output 
listing  appears  in  Table  VII.  Note  the  results  of  mode  changes  and  the  use  of 
command  4  to  generate  special  lists. 


TABLE  VI.  SAMPLE  INPUT  FOR  LISTING  PROGRAM. 


INDEX 

ALL 

GEOMETRY 
ARRAY 
RUN  96 

RUN  122  ID  3 
RUN  122  10  5 
GEOMETRY  IJ  5 
ENO 


3.1.2  Interpolation  Functions 

Two  interpolation  function  subroutines  are  provided  to  aid  the  analyst.  The 
function  ANTRP  interpolates  data  associated  with  endwall  instrumentation,  and 
BNTRP  interpolates  data  associated  with  the  traversing  exit  probe.  ANTRP  and 
BNTRP  interact  with  the  data  base,  retrieving  data  as  needed,  and  both  use 
subroutine  GNTRP  to  obtain  an  interpolated  value.  Subroutine  GNTRP  is  a  user- 
transparent  routine  that  performs  a  weighted  least-squares  fit  of  the  data  and 
returns  the  single  value  for  given  coordinates. 

Subroutines  FCTNL,  CHOLEP,  and  CBSBMP  are  matrix  routines  used  by  GNTRP  to 
determine  the  curve-fit  coefficients.  They  are  also  user-transparent. 
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TABLE  VII.  SAMPLE  OUTPUT  FOR  LISTING  PROGRAM 


«"J 

PSJ 

r< 

o 

2 

O 

ij 

O 

o  2 

•• 

O  2 

V  M 

o  z 

2  ~ 

X 

2  -u 

\ 

2  — 

s 

Q 

QC 

—  a 

QC 

3 

Q  < 

X 

o  < 

X 

o  < 

T 

<  UJUJ 

\ 

<  UJUJ 

V 

<  LUuJ 

\ 

UJUJQCQC 

23 

UJUJQCQC 

z> 

LUlJQCQC 

X) 

« 

« 

QC  ^3 

^ac 

U.ICJ 

*”  X 

ulU  J 

*-u. 

* 

•» 

2»-o 

03  UJ 

^p-w 

duj 

OP-O 

CQLU 

ft 

JP-<\ 

03 

-Jp-<N 

<x> 

JH<\ 

CD 

-J<QC»— 

UJ 

UJ  •  X 

-j<ac*- 

UJ  UJ 1  X 

-J<QCP- 

UJ 

;«J «  r 

<  QCUJ 

QC 

OC  X> 

<  QC  UJ 

QC  QC  13 

<QCUJ 

QC 

^  13 

0CUJO.UJ 

X) 

Z3XZ 

ZUfi.UJ 

3  2>XZ 

XUO.U 

X) 

^X2 

a&£Z 

♦— 

»-2> 

QO.ZZ 

P-  >-23 

oazz 

p-X) 

ZXu< 

< 

<  JZ 

ZIuMt 

<  <uZ 

ZXu4 

< 

<-JZ 

UJUJ^> 

^uu.a 

UJUJP-> 
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TABLE  VII.  (OONT) 
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Function  1 — Endwall  Interpolation  (ANTRP) 

ANTRP  Is  used  to  Interpolate  data  associated  with  fixed  endwall  instrumenta¬ 
tion*  The  formal  input  parameters  for  ANTRP  are  XIN,  YIN,  NRUN,  IDX,  and 
IDAT,  as  defined  in  Table  VIII. 


TABLE  VIII.  ANTRP  INPUT  PARAMETERS 


ANTRP  Parameter 


Definition 


XIN 

YIN 

NRUN 

IDX 


IDAT 


Axial  coordinate 
Tangential  coordinate 
Run  number  (data  set) 
Geometry  I.D.  of 
appropriate  instrumentation 
coordinates 
Data  array  I.D. 
corresponding  to  the 
appropriate  instrumentation 


The  function  software  will  accept  IDX  values  of  5  or  less,  corresponding  to 
endwall  instrumentation  coordinates.  ANTRP  accepts  IDAT  values  of  32,  74,  75, 
or  78  only,  which  correspond  to  endwall  pressures,  hot  surface  temperatures, 
cold  surface  temperatures,  and  adiabatic  wall  temperatures,  respectively. 

ANTRP  also  checks  the  coordinates  (XIN,  YIN)  to  verify  that  they  are  within 
the  valid  endwall  domain.  The  function  returns  a  single  interpolated  value  of 
the  variable  in  C  when  the  statement  C  -  ANTRP  (XIN, YIN, NRUN, IDX, IDAT)  is 
encountered . 

Function  2 — Exit  Plane  Interpolation  (BNTRP) 

The  function  BNTRP  is  used  to  interpolate  raw  and  reduced  data  associated  with 
the  traversing  exit  probe.  The  function  inputs  are  ZIN,  YIN,  NRUN,  IDAT,  as 
defined  in  Table  IX. 


TABLE  IX.  BNTRP  INPUT  PARAMETERS 


BNTRP  Parameter 


Definition 


ZIN 

YIN 

NRUN 

IDAT 


Radial  coordinate 
Tangential  coordinate 
Run  number  (data  set) 

Exit  data  array  I.D. 

(array  dimensioned  (25,  9)) 
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The  function  software  accepts  LDAT  values  of  90  through  107  only,  correspond¬ 
ing  to  exit  data  arrays*  BNTRP  also  checks  the  coordinates  (2IN,  YIN)  to  ver¬ 
ify  that  they  are  within  the  exit  plane  traversed  for  the  run  NRUN.  The  probe 
position  coordinates  are  retrieved  automatically  by  BNTRP.  BNTRP  returns  a 
single  interpolated  value  of  the  variable  in  C  when  the  statement 
C  -  BNTRP  (ZIN, YIN, NRUN, IDAT)  is  encountered. 

The  listing  of  a  program  that  demonstrates  the  use  of  ANTRP  and  BNTRP  is  given 
in  Table  X.  The  program  is  loaded  with  the  data  base  access  subroutines  and 
interpolation  subroutines,  and  reads  cards  specifying  an  interpolation  func¬ 
tion  and  input  parameters  until  an  END  card  is  encountered.  The  program  reads 
input  parameters  from  the  cards,  synthesizes  the  correct  call  statement  of 
ANTRP  or  BNTRP,  and  prints  the  resulting  interpolated  value  for  each  data 
card.  A  sample  input  data  set  is  shown  in  Table  XI  and  the  resulting  output 
appears  in  Table  XII.  This  program  also  uses  the  access  subroutines  DBDMNS, 
DBRUN,  DBSET,  and  DBFCHC  to  retrieve  the  ANTRP  input  data  and  coordinate  ar¬ 
rays,  and  displays  them  for  convenience. 

3.1.3  File  Access  Subroutines 


The  access  subroutines  are  the  building  blocks  of  FORTRAN  programs  interacting 
with  the  data  base.  They  read  and  write  sequential  magnetic  tape  files  and 
create  and  access  a  random-access  data  bank  disk  file.  Several  routines  in¬ 
teract  with  the  random-access  data  bank,  performing  the  specialized  tasks  re¬ 
quired  for  item  storage  and  retrieval.  Table  XIII  summarizes  tasks  performed 
by  the  access  subroutines.  A  description  and  usage  instructions  are  given 
below  for  each  subroutine,  followed  by  comments  on  the  general  structure  for 
data  base  programs  in  Section  3.1.4. 

Subroutine  DBLOAD 


The  calling  convention  for  using  DBLOAD  is  simply  "CALL  DBLOAD".  DBLOAD  is 
called  to  read  the  data  base  from  a  sequential  tape  file  and  generate  a  ran¬ 
dom-access  data  bank  disk  file.  The  subroutine  begins  with  a  call  to  DBSTRT, 
which  initializes  the  master  index  for  the  random-access  file  on  unit  29  and 
writes  the  first  record  of  the  file.  Records  are  written  to  and  read  from  the 
random-access  data  bank  using  the  general  routines  DBPUT  and  DBGET,  respec¬ 
tively.  The  subroutines  DBSTRT,  DBPUT,  and  DBGET  are  user-transparent.  After 
the  data  bank  is  initialized,  records  are  read  from  the  sequential  tape  file 
using  DBREAD,  another  user-transparent  routine.  The  records  are  then  stored 
in  the  random-access  data  bank,  where  a  global  directory  and  set  directories 
map  the  location  of  each  datum  for  future  reference.  Upon  completion  of  the 
loading,  DBDONE  is  called  within  DBLOAD  to  close  the  data  bank  and  write  its 
master  index. 

The  subroutine  DBLOAD  expects  to  find  a  sequential  tape  file  of  the  data  base 
on  unit  27.  TAPE27  and  TAPE29  must  be  specified  in  the  program  statement. 

The  output  resulting  from  a  call  to  DBLOAD  is  shown  in  Table  XIV.  The  number 
of  data  sets  and  the  number  of  variables  stored  for  each  data  set  are  listed. 
Upon  successful  completion  of  loading,  the  final  message  is  printed. 
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TABLE  X.  INTERPOLATION  SOFTWARE  DEMONSTRATION  PROGRAM  LISTING 


MA  INLINE  TO  DEMONSTRATE  USE  OF  FUNCTIONS  ANTRP  AN  J  3nT'P 


10 


LOGICAL  SCAN 

01  HE  NS  ION  CARD (20)  ,  DATA! 5)  ,XCOORD<  300) ,XC 
DATA  NF/5/ 

CALL  DFL2AD 
CALL  OFOPEN 

READ(NR.2,END*97)  CARO 
FORMAT ( 20  A4 ) 

IF (SCAN(CARC  ANTRPt* ) )  GO  TO  20 
IFISCAMCARO  ,'BNTRPi*  )  )  GO  TO  3 C 
IF  1  SCAMCARO  ,  *END$  '  ) )  GJ  TO  09 
IF ( SCAN (CARD,'  ♦'))  GO  TO  1 
HR  ITS (6,10) 

FORMAT (  '  COMMAND  ERROR') 

GO  TO  ] 


A  T  A  (  150 


ANTRP  r  AT  A 
2 G  L*  1 

DO  25  1*1,5 

DATA!  I  )*FGTNUM(C ARD.L, 33 ,L ,NF) 

IF(NF.EQ.O)  GO  TO  97 

25  CONTINUE 
XIN*DATA ( 1 ) 

YIN*DATA(2) 

IDX*OATA( 3) 

NRUN  * DAT A  ( <i  ) 

IDAT  A*CAT  A ( 5 ) 

CALL  DE RUN ( NRUN  , 2 ) 

CALL  OeOMNS( IDATA.MP.MS ) 

CALL  OEFCHC( IDATA.MP.MS, XOATA) 

CALL  DE SET ( 'GEOMETRY' , 2 ) 

CALL  D6FCHCI IDX, 2 .MP.XCOORD) 

WR ITE(6 ,21 ) 

21  FORMAT ( 1H0 . 1 6H ANTRP  INPUT  DATA ,// , 5X 1HX ,5X 1 HY , 8XAH D A TA ) 

DO  22  1*1, MP 

K*2*I-1 

KP1*K*1 

22  WRITE(6,24>  XCODRD ( K ). XCOORD UP  1 >, XDATA ( I ) 

24  FORMAT (2F7.2 ,F12 .2) 

A- AN TRP(X IN .YIN, NRUN, I  OX , I  DATA) 

WRITE (6, 26)  XIN, Y IN, NRUN,  1  OX  ,  I  DATA, A 

26  FORMAT ( / , '  ANTRP  INPUT  X  IN  ,  Y  IN  ,NRUN , I DX , I  DA TA = » , 2F 7 . 2, 31 5 , / , 

1  '  OUTPUT  VALUE*' .F10.3) 
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GO  TO  1 

BNTRP  OATA 
L*  1 

DO  35  1*1.4 

DATA (I  )*FGTNUM(CAR0,L,80,L,NF) 
IF(NF.EQ.O)  GO  TO  97 
CONTINUE 


Z I N*DAT  A ( 1 ) 
YiN»0ATA(2) 
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97 

98 

99 


DATA ( 2 ] 

NRUN*DATA(3) 

IDATA*0ATA(4) 

B*BNTRP(ZIN,YIN,NRUN,IOATA) 

WRITE (6, 36)  ZIN, Y IN .NRUN  ,  I  DATA  ,B 

FORMAT ( '  BNTRP  INPUT  Z IN , Y IN ,NRUN , I  DATA* • , 2F 7 . 2 ,2 1  5 , /, 
1  •  OUTPUT  VALUE*' , FID, 3) 

GO  TO  ,1 
WRITE (6,98) 

FORMAT ( '  MISSING  INPUT  PARAMETER') 

GO  TO  1 
CALL  DEOONE 
CALL  EXIT 
END 
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TABLE  XI.  SAMPLE  INPUT  FOR  INTERPOLATION  DEMONSTRATION  PROGRAM 
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■** 
C  L. 

32 

AN 

TRP 

5 

.5  2. 

,5933 

1 

12 

2 

74 

AN 

TRP 

5 

•  5  2  i 

,5933 

4 

12 

c 

73 

AN 

TRP 

c 

.5  Z, 

.5933 

5 

12 

2 

32 

EN 

TRP 

1 

.2  00 

12.0 

1 

22 

9 

0 

ENC 

Subroutine  DBOPEN 


DBOPEN  opens  an  existing  data  bank  on  unit  29  by  reading  its  first  global 
directory  record  and  reading  its  master  index  into  core  memory.  Once  the  data 
bank  has  been  opened  using  DBOPEN,  data  may  be  accessed  using  the  subroutines 
DBSET,  DBRUN,  DBSTDM,  DBDMNS ,  DBF'CHC,  DBFCHP,  DBSTRC,  and  DBSTRP.  The  calling 
convention  is  simply  "CALL  DBOPEN."  No  printed  output  is  generated  by  this 
routine. 

Subroutine  DBSET  (NAME ,10 PT) 

Subroutine  DBSET  is  used  to  access,  create,  or  delete  a  data  set  in  the  ran¬ 
dom-access  data  bank.  The  subroutine  calling  convention  is  "CALL  DBSET  (NAME, 
IOPT)."  NAME  is  the  eight-character  data  set  identifier,  such  as  GEOMETRY  or 
RUNXXXXX.  IOPT  is  an  integer  corresponding  to  one  of  the  desired  action  op¬ 
tions  outlined  in  Table  XV. 

Option  1  is  useful  for  creating  new  data  sets.  Option  2  is  used  for  accessing 
data  sets  only.  Option  3  is  used  to  ensure  a  unique  new  data  set  is  being 
formed.  For  example,  option  3  is  used  when  the  data  bank  is  loaded  using 
DBLOAD.  Option  4  is  used  to  purge  data  sets  no  longer  needed. 

Subroutine  DBSET  generates  no  printed  output  when  it  successfully  performs  the 
user's  desired  task.  If  option  2  is  elected  and  the  data  set  is  not  found, 
the  message  "DATA  SET  NAME  NOT  FOUND"  is  printed.  If  option  3  is  elected  and 
the  data  set  is  found,  the  message  "DATA  SET  NAME  ALREADY  EXISTS"  is  printed. 
If  the  message  "DATA  BANK  IS  INCONSISTENT"  is  printed,  an  error  occurred  when 
the  data  bank  was  loaded.  This  message  indicates  a  software  problem  in  DBLOAD 
and  should  not  appear  during  normal  operation  of  the  data  base. 

Subroutine  DBRUN  (NRUN.I00DE) 

Subroutine  DBRUN  synthesizes  the  data  set  name  RUNXXXXX  and  calls  DBSET  to 
access  the  data  set.  The  calling  convention  is  "CALL  DBRUN  (NRUN , ICODE) " . 

NRUN  is  the  experimental  run  number  associated  with  a  particular  data  set. 
ICODE  is  the  user  option  associated  with  the  DBSET  call  as  outlined  in  Table 
XV. 
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TABLE  XII.  SAMPLE  OUTPUT  TOR  INTERPOLATION  DEMONSTRATION  PROGRAM 
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TABLE  XII.  (CONT) 
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l.  .50 

2.28 

612.30 

3 .42 

2.25 
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TABLE  XIII.  ACCESS  SUBROUTINE  FUNCTIONS 


Task  Subroutine 


o  Load  sequential  tape  file  into  random-  DBLOAD 

access  data  bank  disk  file 

o  Open  existing  random-access  data  bank  DBOPEN 

o  Access/create/delete  a  data  set  DBSET 

o  Access/create/delete  data  set  DBRUN 

"RUNXXXXX" 

o  Reserve  storage  for  an  array  DBSTDM 

o  Retrieve  dimensions  of  a  stored  array  DBDMNS 

o  Store/overstore  an  entire  array  DBSTRC 

o  Retrieve  an  entire  array  DBFCHC 

o  Store/overstore  part  of  an  array  DBSTRP 

o  Retrieve  part  of  an  array  DBPCHP 

o  Close  random-access  data  bank  DBDONE 

o  Write  data  bank  contents  to  sequential  DBSAVE 

tape  file 


As  with  the  subroutine  DBSET,  there  is  no  printed  output  for  normal  execu¬ 
tion  of  subroutine  DBRUN.  User  error  messages  for  a  given  option  are  printed 
as  generated  in  DBSET. 

Subroutine  DBSTDM  (IV,MP,MS) 

DBSTDM  is  used  to  reserve  storage  in  the  random-access  data  bank  for  a  matrix 
dimensioned  MP  x  MS.  A  call  to  DBSTDM  must  be  preceded  by  a  call  to  DBSET  or 
DBRUN,  which  accesses  a  particular  data  set. 

The  calling  convention  for  this  subroutine  is  "CALL  DBSTDM  (IV,MP,MS)".  The 
user-supplied  arguments  of  DBSTDM  are  variable  identification  number  (IV), 
matrix  primary  dimension  (MP) ,  and  matrix  secondary  dimension  (MS).  The  sub¬ 
routine  adds  the  variable  I.D.  to  the  set  directory  of  the  data  set  accessed, 
reserves  storage  records  for  a  data  array  dimensioned  MP  by  MS,  and  initia¬ 
lizes  the  data  array  members  to  zero  in  anticipation  of  partial  storage  of  the 
matrix.  The  subroutine  normally  generates  no  printed  output.  If  the  DBSTDM 
call  is  not  preceded  by  a  DBSET  or  DBRUN  call,  the  message  "NO  SET  CALLED  BE¬ 
FORE  DBSTDM"  is  printed. 

Subroutine  DBDMNS  (IV,MP,MS) 

Subroutine  DBDMNS  retrieves  the  dimensions  of  a  matrix  stored  in  the  random- 
access  data  bank.  A  call  to  DBDMNS  must  be  preceded  by  a  call  to  DBSET  or 
DBRUN. 

The  subroutine  call  statement  is  "CALL  DBDMNS  (IV,MP,MS),"  where  IV  is  the 
user-supplied  variable  I.D.  The  subroutine  returns  the  primary  and  secondary 
dimensions  associated  with  IV  in  MP  and  MS,  respectively.  Printed  output  is 
generated  only  if  an  error  occurs.  If  the  DBDMNS  call  is  not  preceded  by  a 
DBSET  or  DBRUN  call,  the  message  "NO  SET  CALLED  BEFORE  DBDMNS"  is  printed.  If 
the  variable  number  IV  is  not  present  in  the  current  set  directory,  the  mes¬ 
sage  "UNABLE  TO  LOCATE  IV  IN  SET  NAME"  is  printed. 
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TABLE  XIV.  PRINTED  OUTPUT  FOR  DBLOAD 


45  SETS  BEING  LOADED  IN  BANK. 


SET 
— SST- 

RUN  122  HAS  139 
-RUN - i-U— HAS - 1 21- 

VAR 

-VAR 

IABLES 

IABLES 

SET 

RUN 

91 

HAS 

113 

VAR 

IABLES 

SET 

RUN 

96 
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21 

VAR 

IABLES 

SET 

RUN 

93 
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21 

VAR 

IABLES 

SET— RUN - 94— HAS 

-4-3 

VAR 

IABLES 

SET 

GEOMETRY 
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14 

VAR 

IABLES 

SET 

RUN 

95 
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22 

VAR 

IABLES 

SET 

RUN 

9S 
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43 

VAR 

IABLES 

-  SE  T-RUN-  - 

—  99 

-HAS-  - 

113 

VAR 

IABLES 

SET 

RUN 

105 
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95 

VAR 

IABLES 

SET 

RUN 

107 

HAS 

95 

VAR 

IABLES 

SET 

RUN 

103 

HAS 

95 

VAR 

IABLES 

SET 

Run  - 

-109 

HAS 

113 

VAR 

IABLES 

SET 

RUN 

111 

HAS 

113 

VAR 

IABLES 

SET 

RUN 

112 

HAS 

113 

VAR 

UBLtS 

SET 

RUN 

113 

HAS 

113 

VAR 

lABLtS 

— SET- 

-RUN — 

-57- 

-HAS- 

94- 

VAR 

IABLES 

SET 

RUN 

61 

HAS 

94 

VAR 

IABLES 

SET 

RUN 

67 

HAS 

94 

VAR 

IABLES 

SET 

RUN 

114 

HAS 

113 

VAR 

IaBLES 

-SET- 

-RUN- 

423 

HAS 

121 

VAR 

IABLES 

SET 

RUN 

125 

HAS 

10 

VAR 

IABLES 

SET 

RUN 

131 

HAS 

40 

VAR 

IABLES 

SET 

RUN 

132 

HAS 

51 

VAR 

IAcLES 

— SET — SUN - 133--  HAS 

SET  RUN  149  HAS 
SET  RUN  150  HAS 
SET  RUN  165  HAS 


-SI  VARIABLES 
51  VARIABLES 
51  VARIABLES 
51  VARIABLES 


- SET-RUN — 166- HAS - 5  1  VAR  iABLcS 


SET  RUN  lag  has 
SET  RUN  169  HAS 
SET  RUN  170  HAS 

-SE-T-RUN - 171  HAS 

SET  RUN  172  HAS 
SET  RUN  173  HAS 
SET  RUN  174  HAS 


165  HAJ 
169  HAi 


51  VARIABLES 
51  VARIABLES 
51  VARIABLES 
51  VARIABLES 
51  VARIABLES 
51  VARIABLES 
51  VARIABLES 


-SE-T-RUN — 118— HAS  — 51  VARIABLES 


SET  RUN 
SET  RUN 
SET  RUN 
-SE-T — RUN 
SET  RUN 
SET  RUN 
SET  RUN 


124  HAS 

86  HAS 

87  HAS 
-89- HAS 

176  HAS 

177  HAS 
179  HAS 


-08AN  K--fcGA-£) -COMPLETED  — 


121  VARIABLES 
93  VARIABLES 
93  VARIABLES 
-93  VAR  IABLES 
9  VARIABLES 
9  VARIABLES 
9  VARIABLES 


Subroutine  DBSTRC  (IV,MP,MS , ARRAY) 


Subroutine  DBSTRC  Is  used  to  store  or  overstore  in  the  data  bank  the  entire 
matrix  dimensioned  MP  by  MS.  The  statement  "CALL  DBSTRC  ( IV, MP, MS, ARRAY) " 
must  supply  the  arguments  IV,  MP,  MS,  and  ARRAY.  IV  is  the  variable  identifi- 
cation  number  (I.D.)  to  be  found  in  the  set  directory  or  added  to  the  set  di¬ 
rectory,  and  MP  and  MS  are  the  primary  and  secondary  dimensions  of  the  data 
matrix,  respectively.  The  data  are  stored  in  the  vector  ARRAY.  ARRAY  must 
contain  MP  x  MS  items  or  the  storage  will  not  be  complete.  The  data  should  be 
arranged  in  ARRAY  with  the  primary  index  varying  fastest.  In  other  words, 
elements  of  a  matrix  to  be  stored  in  the  data  bank  under  IV  sized  3  by  2 
should  be  ordered  (1,1),  (2,1),  (3,1),  (1,2),  (2,2),  (3,2)  in  ARRAY  prior  to 
the  DBSTRC  call. 
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TABLE  XV.  DBSET/DBRUN  OPTIONS 


IOPT  Action 

1  Use  data  set  If  found,  create  if  not  found 

2  Use  data  set  if  found,  error  if  not  found 

3  Error  if  data  set  found,  create  if  not  found 

4  Delete  data  set  if  found,  return  if  not  found 


No  printed  output  is  generated  by  subroutine  DBSTRC  normally.  If  the  call  was 
not  preceded  by  a  DBSET  or  DBRUN  call,  an  error  message  is  printed. 

Subroutine  DBFCHC  ( IV,MP,MS , ARRAY) 

Subroutine  DBFCHC  retrieves  an  entire  matrix  from  the  data  bank.  The  user- 
supplied  arguments  are  IV,  MP,  and  MS,  which  are  the  variable  I.D.  and  matrix 
primary  and  secondary  dimensions,  respectively.  The  subroutine  returns  the 
data  in  ARRAY,  with  elements  ordered  with  the  primary  index  varying  fastest. 

For  a  normal  fetch  of  a  complete  matrix,  no  print  is  generated.  If  the  user- 
supplied  primary  and  secondary  dimensions  do  not  match  the  dimensions  stored 
in  the  data  bank  for  IV,  the  message  "DIMENSIONS  DON'T  MATCH"  is  printed.  If 
the  variable  I.D.  (IV)  is  not  present  in  the  current  set  directory,  "UNABLE  TO 
LOCATE  IV  IN  SET  NAME"  is  printed.  If  the  DBFCHC  call  was  not  preceded  by  a 
DBSET  or  DBRUN  call,  "NO  SET  CALLED  BEFORE  DBFCHC"  is  printed. 

Subroutine  DBSTRP  ( IV, IBP , IEP , IBS ,IES . ARRAY, IPI) 

Subroutine  DBSTRP  is  used  to  store  or  overstore  a  part  of  a  data  matrix  in  the 
data  bank.  The  calling  sequence  arguments  IV,  IBP,  IEP,  IBS,  IES,  ARRAY,  and 
IPI  are  supplied  by  the  user.  IV  is  the  variable  identification  number.  IBP 
and  IEP  are  the  beginning  and  ending  primary  indices  for  the  data  to  be 
stored.  IBS  and  IES  are  the  beginning  and  ending  secondary  indices.  The  data 
must  be  present  in  core  in  the  vector  ARRAY.  The  primary  index  IPI  of  the 
dimensioned  vector  AKRAY  must  also  be  supplied. 

DBSTRP  generates  no  print  for  a  successful  partial  store  of  a  vector.  If  the 
variable  IV  has  not  been  initialized  by  a  DBSTRC  or  DBSTDM  call,  the  message 
“UNABLE  TO  LOCATE  IV  IN  SET  NAME"  is  printed.  If  the  DBSTRP  call  was  not  pre¬ 
ceded  by  a  call  to  DBSET  or  DBRUN,  the  message  "NO  SET  CALLED  BEFORE  DBSTRP" 
is  printed. 

Subroutine  DBFCHP  (IV, IBP, IEP, IBS, IES .ARRAY, IPI) 

Subroutine  DBFCHP  is  used  to  retrieve  a  part  of  a  matrix  stored  in  the  data 
bank.  The  arguments  IV,  IBP,  IEP,  IBS,  IES,  ARRAY,  and  IPI  are  defined  as  for 
subroutine  DBFCHP  described  above.  The  part  of  the  matrix  requested  is  re¬ 
turned  in  ARRAY,  which  has  the  primary  dimension  IPI. 
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Subroutine  DBFCHP  generates  no  print  for  normal  execution.  If  the  variable  IV 
has  not  been  previously  stored,  the  message  "UNABLE  TO  LOCATE  IV  IN  SET  NAME" 
is  printed.  If  the  DBFCHP  call  was  not  preceded  by  a  DBSET  or  DBRUN  call,  an 
error  message  is  printed. 

Subroutine  DBDONE 


Subroutine  DBDONE  is  used  to  close  the  random-access  data  bank  file  on  unit 
29.  The  routine  writes  the  total  number  of  records  in  the  bank  on  the  first 
record  of  the  file  and  copies  the  master  index  onto  the  end  of  the  file. 

Subroutine  DBSAVE 


Subroutine  DBSAVE  is  called  when  the  random-access  data  bank  has  been  edited 
and  a  new  sequential  tape  file  of  the  data  base  is  to  be  created.  The  subrou¬ 
tine  reads  each  record  of  each  data  set,  one  data  set  at  a  time,  and  writes 
the  data  to  unit  28  using  DBSTDR.  DBSTOR  is  a  user-transparent  routine  that 
packs  the  data  into  80-character  records  and  writes  the  records  to  unit  28  as 
they  are  filled. 

A  call  to  DBSAVE  requires  no  input  information.  The  subroutine  assumes  that  a 
consistent  data  bank  exists  on  unit  29  and  that  it  may  write  to  unit  28.  Both 
TAPE28  and  TAPE29  must  be  present  in  the  program  statement.  The  printed 
output  generated  by  DBSAVE  is  similar  to  that  of  DBLOAD  and  is  shown  in  Table 
XVI.  The  total  number  of  data  sets  in  the  data  bank  is  printed,  and  each  set 
stored  is  listed  with  the  number  of  variables  in  it  specified.  I 


3.1.4  User-Written  Analysis/Report  Programs 

All  FORTRAN  programs  interacting  with  the  data  base  follow  a  simple  format, 
which  can  be  seen  by  referring  to  Table  XVII  or  a  listing  of  the  maintenance 
program.  First,  the  program  statement  should  designate  all  file  units  used  by 
the  program.  Then  the  data  base  must  be  loaded  into  a  random-access  disk  file 
by  calling  DBIOAD.  If  the  data  bank  exists  as  a  random-access  disk  file  des¬ 
ignated  unit  29,  the  DBLOAD  call  may  be  omitted.  The  random-access  file  is 
opened  by  calling  DBOPEN,  which  retrieves  the  first  global  directory  record. 

Once  the  bank  has  been  opened,  it  may  be  accessed  by  means  of  the  subroutines 
described  in  Section  3.1.3  and  the  interpolation  functions  described  in  Sec¬ 
tion  3.1.2.  To  enable  the  print  messages  associated  with  the  access  subrou¬ 
tines,  the  print  flag  IERPRT  must  be  set  to  1  in  the  FORTRAN  code  and  the  fol¬ 
lowing  common  block  must  be  present  in  the  program: 

COMMON /BUFR/A( 20) ,  LREC,  LA,  IB,  IAS,  IERR0R,  IERPRT. 

A  program  Interacting  with  the  data  base  should  always  be  terminated  with  a 
call  to  DBDONE  to  preserve  the  master  index  for  the  disk  file.  If  a  new  seq¬ 
uential  file  of  the  data  base  needs  to  be  written,  the  DBDONE  call  should  be 
followed  by  a  DBSAVE  call. 
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3.2  GENERAL  OPERATING  PROCEDURES 

The  access  subroutines,  interpolation  software,  and  utility  subroutines  SCAN, 
OOMPAR,  FGTNUM,  and  GETNUM  should  all  be  loaded  into  a  permanent  library. 
Sample  control  cards  used  to  generate  a  library  named  BANK  are  shown  in  Table 
XVIII.  The  last  subroutine  in  the  source  file  should  be  DBDONE.  If  the 
source  code  is  stored  as  a  disk  file,  named  FILE,  the  control  cards  should  be 
preceded  by  an  appropriate  control  card  to  access  the  file,  and  the  FIN  con¬ 
trol  card  should  be  changed  to  read  FTN,I*FILE,R»2.  In  such  a  case  the  source 
deck  is  omitted. 


TABLE  XVII.  GENERAL  FORMAT  FOR  DATA  BASE  PROGRAMS 


Col.  7 

PROGRAM  NAME (INPUT, OUTPUT, TAPE5*INPUT,TAPE6“OUTPUT,TAPE27 .TAPE28 .TAPE29) 

COMMON/ BUFR/A(40 ) , LREC , IA , IB , IAS , I ERROR , IERPRT 

CALL  OBLOAD 

CALL  DBOPEN 

IERPRT  -  1 

FORTRAN  code 

CALL  OBDONE 
CALL  OBSAVE 
END 

TABLE  XVIII.  CONTROL  CARDS  FOR  CREATING  A  LIBRARY  OF  BANK  ROUTINES 
[JOBCARD] 

FTN, R=2 . 

REWIND,  LGO. 

REQUEST, BANK, *PF. 

EDITLIB. 

CATALOG , BANK, RP-999 . 

7/8/9  [eor] 

Source  deck 

7/8/9  [eor] 

LIBRARY(BANK.NEW) 

ADD  (*+DBOONE,LGO) 

FINISH. 

ENDRUN. 

6/7/S/9  (eof ] 
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The  general  control  card  foraat  for  running  MAINT,  BLIST,  or  other  data  base 
FORTRAN  programs  is  shown  in  Table  XIX.  Again,  if  the  program  source  code  is 
stored  as  a  disk  file  called  SFILE,  an  appropriate  control  card  must  precede 
the  FTN  card  and  the  FTN  card  should  be  written  FTN,I*SFILE,R“2.  The  source 
deck  would  then  be  omitted. 


TABLE  XIX.  CONTROL  CARDS  FOR  CDMPILATION/EXECUTION  OF  DATA  BASE  PROGRAMS 
[ JOBCARD  (INCLUDES  TAPE  SPECIFICATION)] 

FTN,R. 

ATTACH, BANK. 

LIBRARY, BANK. 

VSN , TAPE 2 7*nnn , TAPE 28*mmm. 

RE QUEST, TAPE27, PE, NORING. 

RE  QUEST , TAPE  28, PE , RING. 

LGO. 

7/8/9  [eor ] 


Source  deck 


7/8/9  [eor] 


Command s /dat a 


6/7/8/9  [eof ] 


The  program  is  loaded  with  the  data  base  library  BANK,  and  serial  numbers  for 
tapes  27  and  28  are  given.  Unit  27  is  specified  NORING  to  guard  against  writ¬ 
ing  over  the  input  data  base.  The  source  deck  is  followed  by  a  sequence  of 
commands  for  the  programs  MAINT  and  BLIST  or  input  data  for  other  programs. 

To  Implement  the  data  base  system  on  an  IBM  370  computer,  several  software 
changes  must  be  made.  First,  the  file  access  subroutines  and  utility  subrou¬ 
tines  listed  in  the  Appendix  must  replace  the  CDC  subroutines  with  the  same 
names.  Also,  the  error  messages  printed  by  the  subroutines  DBDMNS,  DBPCHC, 


DBSTRP,  and  DBFCHP  must  be  changed  to  print  the  data  set  name  correctly.  The 
CDC  versions  of  the  routines  store  the  entire  set  name  as  NA(1).  The  IBM  ver¬ 
sions  store  the  set  name  in  A(l)  and  A(2).  Thus*  write  statements  that  dis¬ 
play  the  set  name  must  write  A(l)  and  A(2),  rather  than  NA(1).  The  formats 
should  correspondingly  be  changed  from  A10  to  2A4  to  display  the  set  name. 

The  matrix  routines  used  by  the  interpolation  software  contain  DOUBLE  PRECISION 
statements,  which  must  be  changed  to  REAL  *  8  for  equivalent  IBM  function. 

Finally,  the  main  programs  MAINT  and  BLIST  must  be  changed.  The  program 
statements  must  be  removed,  and  the  IF  EOF  statements  must  be  replaced  with 
the  IBM  option  END  ■  NNN  for  read  statements- 


APPENDIX 

REPLACEMENT  DATA  BANK  ROUTINES  FOR 
IBM  SYSTEM  OPERATION 
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phegiduo  J*oi  aum-wT  n 


non 


SUBROUTINE  OBLOAD 

LOAD  THE  BANK  FROM  SEQUENTIAL  FILE 
DIMENSION  B ( 4 0  )  »NA ( 40 ) »C ( 40 ) f NC I  40 ) , NAME ( 2 ) ,0 (  1000 ) 

- eOMMON-ZBUF-R  A  -A-t-  40)  ,LREC  fl  A.  IlbrI  AS  »  I  ERROR-,  IERPRT - 

EQUIVALENCE  ( A ( 1 )  ,  B (  1 )  ,N A ( 1  ) )  ,  ( C ( 1 ) ,  NC ( 1 ) ) 

CALL  D8STRT 

FIND  OUT  HOW  MANY  SETS  AFE  IN  THE  3ANK 

- CALL  DBRf  AUdjNSETSy- - 

WR IT E  ( 6 {40 )  NSETS 

40  FORMAT!  *1*  ,15  ,*  SETS  RE1NG  LOADED  IN  DANK') 

LOOP  THROUGH  THE  SETS 

- IF  INSETS. £0.0)  GO  TO  550  - -  . .  . 

DO  500  ISE T= 1 ,  NSETS 
CALL  DBREA0I2 .NAME) 

CALL  DSREAOd  ,NVAR) 

- WRITE<6,50>  NAME.NVAc  -  -  - - 

50  FORMAT  (*  SET  ',2  A4  ,  ’  FAS*,15,*  VAFlAOLESM 
CALL  DBSET(NAME.l) 

DO  400  I V A R *  1  ,NVAR 

- CALL  DBREAD! 1 , !0N0> 

CALL  OBREAO(I.IPI) 

CALL  DSREAOd, IS!) 

ISIZE*IPI*ISI 

- -|p-^K«ETG-T-i-ieOOl--€0-TO-  150  -  - - —  —  - 

CALL  DBREAOt ISIZE{D) 

CALL  DBSTRC ( I3NQ , IP  I »ISI »D) 

GO  TO  400 

1 50  £Atjj-  09STDHMT1N0 .  IRI  ,  I S I  ) 


CALL 
- CALL 


DO  200  JSI»1,ISI 
IF  ( IPI.GT.1300)  GO  TO  170 
IREAO(IPI.O) 

__  J5TRFd9N0-,l,IPI  ,JSI,  JSI.D.1000T--- 
GO  TO  200 
170  DO  190  JPI-1.IPI.1000 
ILIM*MINO( JPl*999,IPI) 

- CALL  O&REAOdbf  Md-JPI  .O  >  - - - 

CALL  DBSTRP ( IDNQ.JPI  ,ILIM»JSI  .JSI.D.1C00) 
190  CONTINUE 
200  CONTINUE 

AOO  -CONTINUE — - - - - - 

500  CONTINUE 
550  CONTINUE 

CLOSE  OUT  BUFFER 

- CALL-OBDONE -  -- - - 

WR ITE (6j 600 ) 

FORMAT  dOBANK  LOAD  COMPLETED*) 

RETURN 

€*e - 


600 


-StlOR 


CA^T^BS  TRT'TO 
DIMENSION  B (  40)  \  tv  . 

-°WM?ALCNCCR(a{ 1 )*D ( I^k^a 1 1 y  J16  >IAS,IERR0R ,IERPRT 

dIfInE  FILE  2 9( 15600,  *0|U,  IDUM ) 

IERRQR-0 

IERPRT-1 

- BO  10  -I»H40 - 

10  na(!T-o 

LREC-l 

NA<  38>»LREC 

- CALL  DBPUT ( 1 >  A ) - 

IA*0 


NE  DBSTRT 


INITIALIZE 
»NA (  40) 


A  NEW  BANK 
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IB-0 

RETURN 


C 

C 

C 


-tNtr- 


-SUftftflUTINE-  OfrQP^N - - - 

CALL  DBOPEN  TO  OPEN  AN  EXISTING  BANK 
DIMENSION  B (  AO) ,NA (  AC) 

COMMON  /3UFR/  A(  AO )  ,LREC  ,  I  A  ,  IB 


C 

C 

c 

c 

c 

c 

c 

c 

e 

c 

c 


COtt IYALENCE— I  ATI  h»B  M  )  »NA<  1  ) ) 
DEFINE  file  29(15000,  AO.U.ICUM) 
1ERR0RS0 
JERPRT-1 
IA-O 
IB  *0 
IAS-0 

CALL  DBGET ( 1  ,  A) 

~t*EC~NA  t-W - - 

RETURN 

ENO 


IAS. TERROR .IERPRT 


SUBROUTINE  DBSET(NA"E,ICPT) 

CALL  SET  TO  DEFINE  A  DATA  SET 


I  OPT 
1 
2 

—3 - 


VALUES 

USE 

USE 

r  r>  n  n  n 

t  KKy K 


AND  MEANING 
IF  FOUND. 

IF  FOUND. 

_  ~I  F -FOUND-.— 
OELETE  IF  FOUND. 


CFE/.TE  IF  NOT  FOUND. 

E^kOR  if  not  Fnt'N'1, 

— CPF: A-TP-IF-  NOT-  F{jU*fOv- 
RETU^N  if  not  FOUND . 


DIMENSION  NAME ( 2 ) 

DIMENSION  B(  AO)"»NA  (  ACT  - 

COMMON  /BUFR/  A(  AO )  .LREC , I  A, IB 
EQUIVALENCE  ( A ( 1 )  ,5 ( 1 ) ,MA ( 1  ) ) 
IERROR*Q 

“LA*0 - 


I  ERROR  ,  ICRRRT 


IB-0 
ID-1 

10  CALL  DBGET ( ID ,A) 

- IQS*ID - - - 

DO  100  1*1,  31,5 
(NA(I).EQ.O)  GO 


TO 


_ _  _  .150 

NA(I)  .NE.NAMEU  )  )  ^GO^TO^lOO^ ^ 


if 

ioo  Continue 

ID-NA (  AO) 

IF  (ID.NE.O)  GO  TO  10 

- 

LREC-LREC* 1 
NA (  AO) -LREC 
TALL  DBPUT ( ID S , A ) 


90 


TD»tRgC 
00  90  1-1 ,  AO 
NA ( I ) -0 

-  DBPUT ( ID , A ) 


150 


nSu  r 


fm. 

ip.nopT. 


* 

NA 

LR 


IE 


E(J.  2 ) 
“Q.  A ) 


*NAME  (  r 
* ) -NAME  I 
REC*  1 


GO  TO 
GO  TO 


1000 

250 


2) 
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oo  o  rtffinoo 


NA( I*4)»LREC 

-e*tt-e9RUTff  e  t 

I  A*LREC 
I  AS* I A 

NA ( 1 ) -NAME ( 1 ) 


Ah 


40 


200 


TO  30G 


155  DO  160  1*5, 

160  NAU)*0 

NA (  40) *1  A 

- C  Attr-OB  RUT(T*7 AT  - 

GO  TO  250 
IF  (I0PT.EQ.4)  GO 
I A  *NA (1*4) 

IAS* I  A -  - 

CALL  D8GET { I  A , A ) 

IF  (NA(l).NE.NAME(l ) )  GO  TO 
IF  ( NA  ( 2 ) .ME ,NA^P(2 ) )  GO  TO 

TF — rTOPr .  £3v3T~  GO  TOUTO - 

IF  {  IOPT.EO.4 )  Gfj  TO  155 
RETURN 
NA( I )*-l 


105  C 
1C50 


250 

300 


ITE  (6»1C10)_NA*E 


~,2A  4 


hCT-FOL^D-*-)- 


■'rjL^OBPUTXTOTA'r - 

GD*T0  250 

1000  IF  (  IERPRT.NE  .Oi  A 

1  atC-FOfW/rr-T^OO  AT  A-S  E  T 
I ERR0R=1 
RETURN 

1050  IF  (  IERPRT.NE ,0)  WRITE  (6,. 

IOfrO~  FORMATT  *0***OA  TAr  BANK  I  S'  ' 

IERROR-2 
CALL  EXIT 

1100  IF  ( IERPRT.NE.O)  WRITE  (6.1110) 

-mO-FORHAT— (-*il&AT-A-SE-T  -*7  2A4  ,  *  ALREADY  EXISTS  *-f~ 
IERROR-3 
IA-0 
RETURN 
END 


CONSIST  ANT*)- - 


NAME 


SUBROUTINE  D800NE 

CLOSE  A  BANK  BEFORE  END  OF  JOB  OR  OVERLAY 

Thl  ME  NS  I  ON— Bi  40i7NAT-^  Oi - 

COMMON  /BUFR/  A(  40 )  .LREC , I  A  ,  IB , I A S , 1EPRQR , IE RPRT 
EQUIVALENCE  ( A (1 )  ,B ( I )  ,NA ( 1 ) ! 

IERR0R*0 

TA-0 - 

CALL  OBGET (1 ,A) 

NA?  38 )*LREC 
CALL  DBPUT (1 , A) 

-RETURN - 

END 


-e 

C 


»W 


[ON  F ( 4 
[ON-  frC  ■ 


THE  CONTENTS  OF  THE  'BANK 


40  )rNA-f — 40-) — 


MMON  /BUFR/  A(  40 ) ,LREC , I  A , IB , I AS , IERROR , IERPRT 
_  U1VALENCE  (A(i)fNA(l)  •BCD) 

IB-1 
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onnn 


- nrrrx  c — \  onrr  C7i - 

123  FORMAT  MlOUMP  OF  DATA  RANK') 
CALL  DBGET { 1  » A) 

LRECD*NA ( 38 ) 


DO  1000  IR*1 tL^EC 
WRITE  (6,1241  IR 


124  FORMAT  ( *ORECORD* ,15) 

- CAtt-OBSE-H-HHA-1 - - - - 

DO  200  1*1 ,  40,4 
J*I 

IE*I *3 

•JE*^_3 - -  —  -  .  -  -  ■  -  - 

CALL  MVC(F,16,A(  .**  > 

DO  100  K  * 1  ,16 

IF  ( CLCCF ,K , ' A  *  ,  1 ,1 >  4CLC (F ,K  »  •  I  *  ,1,1  )  .LE ,0  * )  GO 

- IF — fCLCf FtKt"*^ *"»  1~»  1  r^CtC  (  F  ,  *R  *  *1,1)  .LE.O.)-  G? 

IF  (CLC(F,K,'S' ,1,1)*CLC(F,K,*Z'  ,1,1). LE.O.)  Gr 
IF  ( CLCi F ,K  ,  '0 •  ,  1 ,1 )  *CLC  ( F ,K  ,  '9  •  ,1 , 1  )  .LE  . 0  • )  GC 
IF  ( CLC( F  ,K ,  '  »  ,1,1  )  .EC.C.)  GC  TO  ICO 

-  5^  CAtL-MVC5  (  FTK  r»-.  *-,  l  ,1>  -  -  - 

100  CONTINUE 

,  WRITE  (6,125)  I,(MA(K),K*I.IE),F,(A(K),<*J,JE) 

125  FORMAT  ( 14 ,4 18 ,3X,4A4 ,3X  ,4G14  .7  ) 

-200 -CONTINUE -  - -  - 

WRITE!  1,124)  IR 
1000  CONTINUE 
_  RETURN  _ 


SUBR 

G 


UTINE  DBGETUREC.BUF) 
NERALIZED  DIRECT  ACCE 


GENERALIZED  DIRECT 
DIMENSION  9UF ( 40 ) 

Re*D-i29^IReC>-9UF - 

RETURN 

END 


ACCESS  DISK  READ  FOR  DATABANK 


SUBROUTINE  OBPUT( IREC,BUF) 

GENERALIZED  DIRECT  ACCESS  DIXK  WRITE  FOR  DATABANK 


29* IREC 


mmmm 


-  DBRUNINRUN, ICODE ) 

jimnw  ji Xt  A  UA  I  A icr  Wflnt  ru^  A  KyW  ”UnD K uN  WN 

IMENSION  NAME(2) 

ATA  NAME  /'RUN  »,*  */ 

L  INCORE (NAME(2) ,4) 


t!name,icooe) 

RETURN 


E 00 EtFH-AL-  PILE— 


DIMENSION  9(40)  ,NA(4C)  ,0 (40  )  ,NC  (  40  ) ,  N  A**E  (  2  )  ,! 
COMMON  /9'JFR/  A(  40  )  ,LREC  .  I  A  .  I S  ,  I A  S ,  1  ERROR  ,  I! 
EQUIVALENCE  ( A ( 1  )  ,9 ( 1 ) ,NA ( 1  )  j , ( C ( 1 ) ,NC(1 > ) 

- RE-W-IN6-2  9 - - - 

CALL  DBOPEN 

FIND  OUT  HOW  MANY  SETS  AFE  IN  THE  SANK 
NSET  S*0 

- - - -  - - - - 

10  00  100  1*1,31,5 
CALL  OSGE  T ( 1 0  ,  C  ) 

IF  (NCm.EQ.O)  GO  TO  15G 
- US-  (Nem.EQ.-l)  GO  TO  100 

SEE  I p  IT  HAS  ANY  VARIABLES 
IC*NC { 1*4) 

CALL  08GET( IC,C) 

— . IF • (NCf5y.E0.01  GO  TO  100  .  *  - 

NSET  S*NSET  S+ 1 
100  CONTINUE 

CALL  09GET ( ID  ,C ) 

.  IF  '(NCf  40)  .E3  .0  )  GO  TO  150 

ID*NC(40) 

GO  TO  10 

150  WRITE(6.160)  NSETS 

160  FORMATHl*  ,15,*  SETS  BEING  STORED  BANK* 

CALL  D8ST0R(1 .NSETS) 

- LOOP- THROUGH-  THE  SETS 

IF  (NSETS . EQ .0 )  GO  TO  55C 

1 0*1 

180  00  500  1*1  ,31  ,5 

- CALL  D96ET f ID  ,C ) 

IF  TnC(I) .EQ.O)  GO  TO  550 
IF  (NC( I) .EQ.-l )  GO  TO  5C0 
NAME}  n*NCj  j  )  _ 

COUNT  THE  VARIABLES 
NVAR*0 
IC*NC(  1+4) 

- ECS  » I C 

220  CALL  OBGET ( IC ,C) 

DO  240  J*5 , 33 , 4 

IF  (NC(J) .EQ.O)  GO  Tp  2r0 

- IF  (Ne(  J^l  >*NCt J+2> .15.0)  GP  TO -240 

NVAR  *NVAR ♦ 1 
240  CONTINUE 
IC*NC(40 ) 

250  IF  (NVAR.EqTO)  GO  TO  500 
WRITE (6,260 )  NAME  »NV AP 

260  FORMAT  (•  SET  ',2 A4,’  HAS*, 15,*  VARI ADLES  * ) 


D ( 1  COO) 
tv  DPT 


CALL  DBSTOR  2  »NAME ) 
CALL  DBSTOR (1, NVAR) 


VW>  '  '  I  W  V  1  «  V  ^  t  M  W  I  V  "  1  I 

^ALyc0BSET(NAME,2> 

300  CALL  DBG6T ( IC ,C ) 

??  ?N?(JK£q*6)  GO  TO  410 
IONO«NC( J) 

— iuzmu) - 

IF  ( !PI*ISI.LE.O>  GO  TO  4G0 


NC  VARIABLES 
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c 


CALL  gBSTQR |I ,I0N3) _ 

CALL  08ST0RU  !lSI) 

4p-X-(4s  I  III?  GfIhKXH-  G  f?  -T  3  31  e - 

CALL  DBFCHCf IDNQ.IPI ,!SI ,0) 

CALL  OBSTQR(  I  $IZE»D  ) _ _ _ 

GO  TO  400 
DO  390  JS I  *  1  *  I  SI 
IF  (IPI.GT.1300)  GO  TO  320 
-CAtL'  -08F-CHf><-f 1  vl P I  » JS I .  JS I ,  C  .  I  CCC  )  - 

CALL  DBSTQR ( IPI ,0) 

GO  TO  390 

00  350  J»I«1  ,IPI ,100C 
ILIM«MIN0{ JPI ♦999, IPI 5 

CALL  D3FC4P(  IDNO.JPI  ,  ILIM  ,  JS I , JSI ,0 , 1 LCO ) 


ISZ«1LIMU-JPI 
CALL  DBSTQRUS 


ISZ.O) 


_  *TfNUE— 

CONTINUE 
CONTINUE 
IC*NC ( 40 ) 

IF  IIC.NEvICSF  GO  T0~300  - 

CONTINUE 

END  OF  VARIABLE  LOOP 
CONTINUE 

CALtrOBOETUO-ve-I - - 

IF  (NC(40) .E3.0)  GO  TO  550 
ID»NC(40) 

GO  TO  180 

- E-NO-  OF  -SE-T-  tOOP - - - - 

CONTINUE 

CLOSE  OUT  BUFFER 
CALL  DBSTQR ( 0  »0 ) 

MRITE ( 6,600-) - - 

FORMAT  ('OBANK  STORAGE  CCMPLETE  C  *  ) 
CALL  DBDONE 
REMIND  23 

RETURN - 

END 


SUBROUTINE  DBSTORtN.O) 

CHINE  DEPENDENT  ROUTINE  FOR  OUTPUT  TO  SEQUENTIAL  FIL 
SION  D(N-)-rC-AR0F20P -  — 

IhM'WWtt'c 

COPY  D  INTO  CARO 


nnn  !  orinn 


IC*0 

- 5e-iP-He>r€ft»-»0  -RETURN- 


GO  TQ  10 
100  IF  ( IC.EQ.Q)  RETURN 
ICP-ICM 

- BQ-MO-f**  .20- - 

110  CARO ( I ) *0 

WRITE ( IOUT  .20  )  CARO 
IC*0 

- RETURN - 

END 


0)  RETURN 


- SUBROUTINE  -OBREAOfN.  01  - 

MACHINE  DEPENDENT  ROUT  INC  FOR  INPUT  FPQi-  SEQUENTIAL  FILE 
DIMENSION  0 ( N ) *  CARO ( 20 ) 

DATA  IC  /20/.IIN  /2  7/ 

- io»o - — - - - 

IF  <  IC.NE.20J  GO  TO  50 
i0  READUIN.20)  CARD 
20  FORMAT  (20A4) 

- ic*0 - - 

50  I D* I D ♦ 1 
IC*IC*1 

D (  ID)*CARD( IC) 

- 1  ^  f-j  0  .  E0-.  -1 )— RE  TURN - -  -  - 

GO  TQ50  *20)  G°  Tn  10 


- FUN&T-IUN— FGTNUMTC4R  DyNS T APT-, NEND  »K  .NCHAR) - 

DIMENSION  CARD! 1 ) 

FGTNUM*0. 

CALL  GETNUM  (CARD, NS  TART,  NEND,  XDUH,  K  ,  NCHAR  ) 

— — IF-  TNCHARrNEvO)-  PGTNU“»XDUM  .  . .  . . 

RETURN 

ENO 


SUBROUTINE  GETNUM (CARO ,N START, NFND ,X NUM.K ,NCHA<? ) 
DIMENSION  CARO(l) 


REAL  MINUS 
DATA  BLANK/' 
1  PLUS2/'* 

-2 - COHNA/', 

NCHAR*0 

XNUM-0.0 

KL*NENOM 

IF (NENO.LT «N 

Nl-NSTART 

N2-NEND-1 

I-NSTART-1 


•/,Z/'Z  '/.DEC/'.  '  /  .PLUS  1  / '  F.  V, 

•/,  MINUS/'-  '/»  E QU A L 1 / ' »  '/,  ECU AL2 / 1 =  V 
-+/-&■  t *r7P/*o v - 


-00~Tfr"3tK 


IF? CLCfCARO  »K »Z  » 1  *  1 1 .GT.0.0 .OR . CLC (C  ARD.K . DEC » 1 » 1 ) . EC .0.0. OR. 
1  CLCTCARO.K.PLUSl.l  .D.EO.O.O.OR.CLC  (CARD  ,K .PLUS2 , 1 , 1 ) .EO .0 . 
2A2l!LAiiCAR6,K,MINUS.l,l).EC.O.C)  GO  TO  100 _ 


85fSWrtf5$2<”^fa  ** 

K»I  +  1 


IF( ( CLC (CARD. I .BLANK ,1,1) .EC .0. 0 .OR .CLC ( CARD. I . EQUAL  1 *1.1) .EC. 


2-.tjRretete*Rthri^QU  At2,  iiit.ts.0.0.  or. etc  <-e  arc  rt  .ccmma ,  i ,  i  j-.l-o 

1  ( CLC (CARD ,K,Z, 1,1) .GT,O.O.CR.CLC(CARQ»<»DEC»1 ,1 ) .EO.O.O.OR. 

2  CLC (CARD, K  .PLUSl.l  ,1) .  EQ  .  C  .0 .  CP. .  CLC  (CARD  ,K  ,PLUS2 , 1 ,  1 )  .EQ.G. 

■SrOR  •  CLC  ( C  ARuyK-fMT-NtlSy  I  ri i~.  ECiCr.O  Hh  G1}-TO-tO© -  - 


50  CONTINUE 
GO  TO  300 
100  CONTINUE 

- me  «~o 


IF(C LC (CARO, K, DEC, 1,1) .EO.C.O)  1DEC  =  1 
IF(K.EO.NEND)  GO  TO  200 


kpi  *  «n 

- 00-150-r*RT»l  »N€NC - -  ~ - - 

KL»I 

IF(CLC(CAR0,I,Z,1,1 ) .GT.C.O.PR. ( CLC ( C  APD , I ,OEC  >1,1) .f  ..0 . 

1  IDEC.EQ. 0) )  GO  TO  125 

- 60- T 0-200 -  -  - - 

125  CONTINUE 

IF(CLC(CAR0,I,DEC,1,1) .cQ.O.C)  ITEC  *  1 
150  CONTINUE 

-  •  KL*KL*1 — . - - - 

200  CONTINUE 

NCHAR  *  MINO<  20.KL-K ) 

00  250  1*1 ,5 

- F  ICLO(  H  BLANK . . - . - . — 

250-  CONTINUE 

CALL  MVC5(FIEL0,21-NCHAR ,CAFC,< , NCHAR) 

CALL  INC0RE(FIELO,2O) 

. REAOt5yi-03OT-XNUM -  -  -  - 

IF(KL.GE.NEND)  GO  TO  3C0 
Nl-KL 
N2*NEN0-1 

- DQ-400-I^t+l-»N2 -  -  - - — -  -  - 

IFCCLC (CAR 0,1, BLANK ,1,1) .NE .C.O)  GC  TO  410 
400  CONTINUE 
GO  TO  300 

4L0"  TF  (CLC  (C:AR-^rI  ,ETL.  T  )-. NE  .0.0  » ANQ  Ztt  TCARD'i  1  ,0,1,1 )  .NE  .0.0) 

X  GO  TO  300 
IP1-I+1 

DO  420  K*IP1 ,NENO 

- IF  ( CLCfCARO ,K .BLANK , If  1 )  ,NE .C .0  )  GO  TO  430 

420  CONTINUE 
GO  TO  300 

430  IF ( CLC (CARO, K  ,Z,  1,1) .LE. 0.0 .AND . CLC ( CARO, K .PLUS  1 ,1,1 ) .NE . 0 .' 

- -X - CtCt«^0-,*CTNLUS?,  l  ,l)  ,t)E  .0.0  ./  ■■■ 

X  GO  TO  300 
KL-K+1 


.ANO.CLCICAPPvK.-lNUS, 1  ,1  )  ,‘Ji 


NENO)  GO  TO  440 


IF 

00'450"I*KPI,NEN0 
KL  ‘ 


(K .EQ. 

ww- 


IF LC (CARD , I ,Z ,1 ,1 ) .LE.C.O)  GO  TO  44C 
-4-50-COfTT-ltWg - 


440 


KL«KL*1 
CONTINUE 

N^HAR1»MIN0120,KL-K ) 


FIELD? 1 1  * 
460  CONTINUE 


INO( 


BLANK 


CALL  MVC5(FIEL0,21-NCHAR1,CARD,KL  -NCMAR 1 ,NCH AR1 ) 
. . 4 --ft Dy2©) - 


CALL  MVCB(FIELO, 

mm 

NCHAR»NCHAR*NCHARIM 
300  K-KL 


61 


RETURN 

1000  FORMAT (F20 .0 ) 

1001  FORMAT ( 120 ) 


LOGICAL  FUNCTION  SCAN(CAH) , LORD ) 

- truGFOAt^l—C  A^O-f-l  >-»  WORD  (■  1 T . ©CL AR -  - 

LOGICAL  COMPAR 
DATA  OOLAR  /•$'/ 

SCAN-. FALSE. 

- 00  10  LEN-1,80 

IF  (CLC(W3R0,LEN«1,DCLAR,1,1).EC.C.>  CC 
10  CONTINUE 
RETURN 

-20  StATt^COMP-A-Ri-C ARO ,  ?OtLE  •■n  WCRf) ,  J ) 

RETURN 

ENO 


TO  20 


LOGICAL  FUNCTION  COf'PAF  (  C  ,  I  ,  L  ,N  ,CU  AF.  ,  J 1 
LOGICAL  *1  C (  1 )  » CHAR (  1  >  ,  ?  LK  ' 

DATA  BLK  /*  •/ 

- NUM»N- -  -  -  - - - 

00  5  K-l.N 

IF(CLC(CHAR,<  ,3L<.1 ,1) .EO.O. )  NLT-NUM-1 
5  CONTINUE 

- IL-L-NUM-+1— . . 

IF(I.GT.IL)  GO  TO  1C 
00  15  K-I.IL 

IF(CLC(C»K»8LK,1,1)  .EQ.O.)  CC  TC  15 

- K1~K - - - - -  - . . 

K2-1 

IL1-IL-1 

20  IF(K2.GT.N)  33  TO  25 

— lF<CLCjCHAR,X2iaLK*1,l).NE.C. )  CO  TO  30 
GO  TO  20 


30  IL1- IL1  +  1 

- 40— EFf  <-l~*-GrT  .-I  tlh -60--T0-1 5 


IF(CLC(C»K1«3LK»1.1)  ,NE.O.)  CC  TO  45 
K 1 -K 1 *1 
GO  TO  40 

-  45- ! Ff  CL  C ( er*-T>-e HA R- , K 2,1  K  NF-. O .  >  GO  TC-  15 

Kl-K 1+1 
K2-K2+1 
GO  TO  20 

-"ira - 

25  J«K1 

COMP AR-. TRUE. 

- RtfiiRN - - - 

10  J-I 

-  COMPAR-. FALSE. 

RETURN 


*U.S.Qo».rnm*nt  Printing  Offlct:  19*1  -  559003/201S 


